Local Empathy provides Global Minimization of Congestion in Communication Networks 
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We present a novel mechanism to avoid congestion in complex networks based on a local knowledge of 
traffic conditions and the ability of routers to self-coordinate their dynamical behavior. In particular, routers 
make use of local information about traffic conditions to either reject or accept information packets from their 
neighbors. We show that when nodes are only aware of their own congestion state they self-organize into a 
hierarchical configuration that delays remarkably the onset of congestion although leading to a sharp first-order 
like congestion transition. We also consider the case when nodes are aware of the congestion state of their 
neighbors. In this case, we show that empathy between nodes is strongly beneficial to the overall performance 
of the system and it is possible to achieve larger values for the critical load together with a smooth, second-order 
like, transition. Finally, we show how local empathy minimize the impact of congestion as much as global 
minimization. Therefore, here we present an outstanding example of how local dynamical rules can optimize 
the system's functioning up to the levels reached using global knowledge. 

PACS numbers: 89.75.Fb, 02.50.Ga, 05.45.-a 



I. INTRODUCTION 

Complex communication networks have recently attracted 
a lot of attention from scientists due to the discovery of the 
topological features of real communication systems such as 
the Internet (2), The structure of these communication sys- 
tems is efficiently described by a graph in which nodes repre- 
sent routers and edges account for the communication chan- 
nels. However, the structure of these graphs is far from being 
purely random. Quite on the contrary, they typically show a 
scale-free (SF) distribution for the number of communication 
channels departing from and arriving to a system's element. 
The use of modern complex network theory J3-01 together 
with tools inherited from non-equilibrium statistical physics 
I!] have allowed to study the dynamical properties of such 
communication systems. In particular, this approach has been 
successfully applied to the study of the structural evolution 

0] of the Internet, its navigability ifsl— TToll and dynamical 
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A lot of the recent literature on communication networks 
has tackled the critical properties of their jamming and con- 
gestion transitions i20M29ll . These studies have focused on the 
design of efficient routing strategies that, on one hand, provide 
with short delivery times and, on the other hand, avoid the on- 
set of the congested state in which the load of packets in the 
system increases, thus causing the failure of information flow. 
It has been shown that finding the best suited strategy depends 
strongly on two main features: the topological patterns of the 
particular network and the load of information on top of it. 
Regarding the first of these two issues, a number of routing 
mechanisms have been studied on different structures |29l432 
allowing to design resilient network backbones 03143511 . 

Many of the routing policies proposed so far rely on the 
(static) structural properties of the communication network. 
Examples of such policies are biased random walks l36l[37ll . 
shortest-path ll38i[39tl and efficient-path |(40f| schemes. These 



routing mechanisms can be conveniently reformulated to in- 
corporate the information about the dynamical state of the 
system, i.e. the congestion state of routers. This allows to 
dynamically change the paths followed by information pack- 
ets in order to bypass those over-congested routes. In this 
line, congestion-aware schemes have significantly improved 
the performance of biased random walks fiUl . shortest-path 
JUIH] and efficient-path JH] routings. In addition to the de- 
sign of efficient routing protocols, several strategies to avoid 
congestion have been implemented. Remarkable examples 
of these strategies are the implementation of incoming flow 
rejection 1451 14611 and packet-dropping mechanisms 114711 for 
avoiding the congestion of single nodes, or the addition of a 
router memory to avoid packets getting trapped between two 
adjacent nodes ll48ll . 

All the above studies have assumed that both network topol- 
ogy and the mechanisms to avoid congestion are static {i.e. 
neither topology nor the routing strategies change). However, 
this approach neglects that, even for the same graph, the opti- 
mal routing po licy de p ends str ongly o n the state of congestion 
of the system Ji2, Hj, |H |H iMl ■ Therefore, in order to 
balance correctly the congestion in a communication system 
it seems appropriate to allow the elements (routers) to switch 
to the best suited strategy to avoid congestion given the in- 
stant traffic conditions. In this article, we propose an adap- 
tive mechanism that allows nodes to choose their individual 
strategies instead of imposing a common policy. In this adap- 
tive protocol routers exploit their local information about the 
congestion state of the system to decide whether to accept in- 
coming packets. First, in section [III we introduce a minimal 
routing model without any adaptive mechanism that allows 
us to unveil the role of rejection when it is externally tuned. 
In section [Til] we will consider that each router can adopt its 
own rejection strategy and make some analytical derivations 
about the optimal strategic configuration to avoid the conges- 
tion onset. In section [IV] we will implement our first adap- 
tive mechanism and show that when nodes are allowed to dy- 
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namically adapt their own strategy, while only being aware of 
their own congestion state (myopic case), the onset of con- 
gestion is shifted to a larger critical load (with respect to the 
static algorithm introduced in section HJi. This improvement 
is due to the self-organization of the strategies of nodes into 
degree-correlated configurations. However, we will show that 
the delay of the onset of congestion comes together with a 
sharp, first-order like, transition that provides no dynamical 
signals about the onset of congestion. Finally, in sectionlVlwe 
show that when nodes are aware of the congestion state of its 
nearest neighbors and empathize with them, it is possible to 
recover the former large critical load together with a smooth 
phase transition, avoiding the uncertain scenario of the my- 
opic adaptive model. More importantly, we will show that 
tuning conveniently the degree of empathy between routers it 
is possible to recover, through a local mechanism, both the 
congestion levels and the rejection patterns provided by the 
global minimization introduced in section Hill 



II. MINIMAL TRAFFIC MODEL 




FIG. 1: (Color online). Phase diagrams, p(p), of the minimal traffic 
model using different values of the rejection rate 77. The inset shows 
the existence of different critical values p c when varying 77 



Let us start by introducing the minimal traffic model in 
which the adaptive algorithm will be implemented below. In 
this model, we consider the transfer of information packets 
between adjacent routers as a probabilistic event. Inspired by 
ll45l , l46ll . we consider a set of stochastic equations for describ- 
ing the time evolution of the queue length of the nodes at some 
time t, q* = {<?*}. The queue length of a given node, q\, can 
either increase or decrease due to several events. First, at each 
time step and with probability p, a new packet is generated 
being added to the queue of the node. Second, at each time 
step each node tries to send a packet in its queue to any of 
its first neighbors. This packet can be rejected by the chosen 
neighbor with some probability 77. If the packet is accepted, 
it may be removed from the system with certain probability 
p. These two latter events mimic the effects, although with 
some important differences, of an active queue control strat- 
egy as the random early detection (RED) 15111 present on Inter- 
net routers and the arrival of the packet to its final destination 
respectively. Following the above ingredients we can write the 
time-discrete Markov chain of the minimal traffic model as: 
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where A^ represents the (i, j) term of the adjacency matrix of 
the network substrate and Q(x) is the Heaviside step function 
(Q(x) = 1 if x > and Q(x) = otherwise). Since our 
network is undirected and unweighted, the adjacency matrix 
is defined as Aij = Aji = 1 if nodes i and j are connected 
and Aij = Aji = otherwise. The quantity ki is the degree 
of a node i — ki), i.e. the number of routers con- 

nected to it. The right-hand-side of equation (Q]i contains two 
terms accounting for the incoming flow of packets that arrive 



to the queue of node i, namely, p (accounting for the external 
load of packets) and the first sum (accounting for the arrival 
of packets from its first neighbors). On the other hand, the 
second sum in equation (Q]i accounts for the probability that a 
packet from i is delivered to a first neighbor. 

The set of equations (Q]i are solved starting from a zero con- 
gestion state: q® = Vi. The evolution of the system is mon- 
itored by means of the following order parameter l23ll : 



P(t) = lim 

T— s-oc 



Q{t + T)-Q{t) 
P T 
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where Q(t) is the sum of all the queue lengths at time step 
t, Q(t) ~ J2i=i 1i- The stationary value, p, of the above 
order parameter is bounded (0 < p < 1) and describes the 
dynamical regime in which the system ends up. Namely, p = 
indicates that the system is able to balance the incoming flow 
of new packets with a successful delivery of the old ones. In 
this case the system is said to operate in the free-flow regime. 
Instead, when p > the above balance is not fulfilled and the 
queues of the nodes increase their size in time at a rate p ■ p. 
In this latter situation the system is in the congested phase. 

We have studied the behavior of the order parameter p by 
taking the rate of packet creation p as the control parameter. 
The arrival-to-destination probability is set to p = 0.2 as the 
usual value found in the Internet JT|] . The corresponding phase 
diagrams are shown in Fig.[T]for several values of the rejection 
probability 77 using a SF network of N = 5000 with P(k) ~ 
k~ 2 2 . As observed in the figure, the transition from free-flow 
to congestion occurs in a smooth way at low values of p being 
the critical point p c = 0.02 for 77 = (no rejection). However, 
as the rejection rate 77 increases the value of p c decreases and 
p increases faster (see inset in Fig.Q]). 
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m. ANALYTICAL APPROXIMATION OF GLOBAL 
CONGESTION MINIMIZATION 

The above results question the convenience of implement- 
ing a rejection mechanism in routing models. However, the 
bad performance of this rejection mechanism relies on the 
homogeneous distribution of the rejection rates across the 
routers of the network. We now explore the general situa- 
tion in which the individual rejection rates are independent. 
Therefore the set of equations (Q]i transforms into: 
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This new set of equations is now used to determine the optimal 
set {rji} so that congestion is minimized for a given value of p. 
To this aim, we first use two assumptions: (i) the nodes have 
reached a stationary state, q\ +1 = q\ Vi, and (ii) the queue 
length of nodes is nonzero, 0(<?f) = 1 Vi. These provisos 
admitted, equations (0 turn into the following set of equations 
for the rejection rates of the routers 
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Now we make use of the annealed approximation of the adja- 
cency matrix i52U54ll : 
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where (k) is the average degree of the network ((k) ~ 4 in our 
case). Introducing the annealed expression (0 into equations 
dUi we obtain: 
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where (fc(l — 77)) = J^. ~ Vj)/N- Equation (0 clearly 
shows that the larger the degree of a router the larger its re- 
jection rate. Therefore, from this expression we observe that 
a non-homogeneous distribution of rejection rates across the 
routers is beneficial to assure the free-flow condition (and thus 
to delay the onset of congestion). We can calculate the expres- 
sion of the rejection rate by computing the value of (k(l— rj)}. 
From equation (0 we obtain: 

(k(l - 77)) = -4- [<*(1 - V)) - P(k)] , (7) 
and finally we have: 

(k(l-ri)) = ^{k). (8) 

Therefore, the rejection rate of a node with connectivity k% 
reads: 
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FIG. 2: (Color online). Rejection rates of nodes as a function of their 
degree r]i(ki) as obtained from equation ((9). The curves correspond 
to different values of the external load of information p. 



As anticipated above, expression (O shows that the rejection 
rates of nodes should depend on their degrees rather than be- 
ing externally set to a constant value. In Fig. |2]we apply equa- 
tion (0 to plot the rejection patterns corresponding to different 
values of the external load p. As shown, 77; decreases with p 
and increases with fcj . 

The assumptions made in order to obtain equation (0 point 
out that its validity, for all the nodes, should be restricted to 
the proximity of the critical point p c . First, for p < p c many 
of the queues are zero (invalidating assumption (ii)) thus mak- 
ing the rejection rate imposed by equation (0 too restrictive 
for the real traffic conditions. On the other hand, for p > p c 
assumption (i) does not hold for all the nodes. This is mani- 
fested by the prediction of negative rejection rates, 77^ < 0, in 
equation (0 for those nodes with low connectivity. In prac- 
tice, the impossibility of displaying negative rejection rates 
fixes their rejection rate to ?/,; = 0. However, those nodes with 
large enough connectivity can still avoid congestion by means 
of positive rejection rates as described in equation (0 (see 
Fig. 0. Following these arguments, we can estimate the exact 
value of p c as the maximum value of p for which 77^ > for all 
the nodes in the network. In particular, given that, for a given 
p, the value of 77^ increases with fej we obtain p c imposing in 
equation (0 that those nodes with the minimum connectiv- 
ity, ki = k m i n , have 77^ = 0. Since in our case k m i n = 2 
and (k) ~ 4 we obtain p c ~ 0.1. Therefore, by externally 
fixing the rejection rate of each node as dictated by equation 
(0 we can assure the permanence in the free-flow phase up to 
Pc^O.l. 



IV. MYOPIC ADAPTABILITY 

The minimal traffic model introduced in section [TT] shows 
that system's performance deteriorates as soon as rejection 
rates are uniformly set in the system. However, in section UTTl 
we have shown that a non-uniform configuration for the re- 
jection rates shifts the critical load to larger values. However, 



4 



this non-uniform configuration has been externally imposed 
and derived analytically following different assumptions. A 
correct derivation of the optimal configuration would imply, 
on one hand, a more sophisticated calculation and, on the 
other hand, a complete knowledge of the architecture of the 
network. This latter condition makes unfeasible the external 
tuning of the individual rejection rates. 

In order to overcome the need of global knowledge about 
the topology of the network we now introduce an adaptive 
scheme based solely on the local information available to 
nodes. In this adaptive setting we will allow nodes to choose 
their own rejection rate so that the dynamical state of a node 
will be described by both q\ and r\\: 



^ = it 
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The individual choice of each instant value 77* aims at oper- 
ating at the optimal regime as given by the external parame- 
ters p and p. To this aim, each node chooses its own rejec- 
tion rate for the following time-step attempting to reach an 
optimal queue length, q opt = p/p, so that traffic is homoge- 
neously distributed across the network. To this end, a node 
raises or decreases its own rejection rate depending on the 
deviation of its instant queue length from the optimal queue, 
A* = q\ — q opt . This rationale mimics a myopic behavior by 
which, regardless of the congestion state of the system, nodes 
are allowed to close the door to new packets while decreasing 
their respective queues. To incorporate this adaptive behavior 
we couple equations ( TTOb with the following evolution equa- 
tions for the set {?7*}: 



= 



1 



1 + cxp f-^A 



(11) 



This evolution rule takes the form of the saturated Fermi func- 
tion so that congested nodes, q. L > q opt , will tend to total rejec- 
tion, rjl +1 — > 1, whereas those under-congested will open the 
door to new packets, r/* +1 — > 0. The velocity of the transition 
from these two regimes is controlled by (3 since it accounts for 
the reactivity of nodes to congestion. Note, that r]* +1 = 0.5 
will be adopted whenever q\ = q opt . 

The adaptive equations ( fTTT i allow for abrupt changes in the 
rejection rates between two consecutive time steps. Thus, we 
also explore a different formulation: 
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in which the rejection rates evolve smoothly. Rule ( TT2T > is 
completed by assuring that 77, remains bounded so that < 
Vi < 1- In the above equation (fT2l . j3 acts as the inverse of 
the time between two consecutive time steps of the adaptive 
dynamics. Therefore, in the continuous time approximation 
of equation (fT2l , the derivative of the rejection rate is equal to 
the difference between the instant queue length and its opti- 



mal value, i.e. A* = q\ — q opt . 
q\ = q opt a router will adopt r/,* +1 



Note that in this setting when 
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FIG. 3: (Color online). (Top) Phase diagram p(p) for the myopic 
routing models A (squares) and B (diamonds) and for the minimal 
routing model (circles). (Bottom) Average rejection rate (77) as a 
function of p of the former three routing schemes. 



In the following we will use the two formulations for the 
myopic adaptive model and show that the results are qualita- 
tively the same. Namely, we will call model A to equations 
( [Tol l and ( fTTT ). and model B to the formulation using equa- 
tions ( [Tol l and ( fT2b . Note that in both models the parameter 
(3 controls the reaction speed of nodes to congestion. In this 
direction, our numerics have shown that by changing j3 one 
basically controls the duration of the transient time before the 
stationary distribution of the rejection rates is reached. In the 
following, we set /3 = 10 and f3 = 10 -2 in models A and B 
respectively. 

In the top panel of Fig. [3] we show the phase diagram, pip), 
of the myopic adaptive model with the two formulations. As 
observed, in both formulations the myopic model displays an 
abrupt, first-order like, transition from the free-flow to the 
congested state. Moreover, in Fig. [3] we have also plotted the 
phase diagram of the minimal model when 77 = 0, i.e. its 
most congestion-resilient version, to show the improvement of 
myopic adaptability by shifting the jamming transition from 
p c = 0.02 to p c ~ 0.1. This value for the critical load is ex- 
actly the same as the one predicted in section [III] by means of 
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FIG. 4: (Color online). Distribution of the individual rejection rates 
rji across degree-classes for several values of p in the myopic routing 
models A (top) and B (bottom). 



organizes homogeneously around r) ~ 0.4. However, when p 
increases the rejection rates of low-degree classes decreases 
while hubs start to close their doors progressively as p in- 
creases. For model B the microscopic configurations adopted 
as p increases are similar regarding the behavior of high- 
degree nodes. However, in this latter scenario low-degree 
nodes remain accepting incoming packets up to the congested 
state. These two figures show that the two different internal 
dynamics (showing different microscopic organizations ) lead 
to the same macroscopic result: the delay of the onset of con- 
gestion. 

Let us highlight that the delay of the congestion onset in this 
myopic adaptive setting again contradicts the results obtained 
for the minimal routing model in which, even a small (ho- 
mogeneously distributed across routers) rejection rate leads 
to an increase of the congestion in the system. Quite on the 
contrary, the myopic adaptive model points out the same idea 
concluded from the global minimization of congestion: a hi- 
erarchical (degree-based) organization of the rejection rates 
by the system is strongly beneficial to avoid the congestion 
of the system. However, from figure [4] it becomes evident 
that the strategies self-adopted in the myopic adaptive set- 
tings are clearly different than the ones obtained in section 
Ulllfrom equation (0 when congestion was minimized using 
global knowledge. Although in equation (0 the value of the 
rejection rate increases with the degree of the node (as in the 
myopic setting), the evolution with p is quite different. Thus, 
although the critical load has been shifted to the same value as 
the one found in section [HI] the self-organized patterns of the 
rejection rates in the myopic settings reveal a clearly different 
scenario. 



V. EMPATHETIC ADAPTABILITY 



the analytical approximation using global knowledge. Thus, 
the myopic adaptive model, equals the delay of the congestion 
onset obtained by minimizing congestion globally. 

To analyze the roots of the resilience of the myopic adaptive 
routing to congestion we have plotted in the bottom panel of 
Fig. [3]the mean value of the rejection rate, (77) = J2iLi Vi- 
lli this case we observe that models A and B display the same 
pattern after the sharp transition to congestion, i.e. the sudden 
closing of all the doors in the network thus causing the abrupt 
transition to p ~ 1 as soon as p > p c . On the other hand, 
the configurations adopted by both models before the onset 
of congestion, p < p c , are quite different: While in model 
B (77} ~ 0, for model A a significant part of the population 
adopts rji > 0. Surprisingly, in this latter setting the average 
rejection rate decreases as we approach the critical point, p c . 

To have a deeper insight about the microscopic configura- 
tions that allow to delay the onset of congestion we show in 
Fig. |4]the set of individual rejection rates of nodes {rji} ranked 
according to their degrees. In both models A and B, the cor- 
relation between rji and fc, is clear since all the routers within 
the same degree-class display similar rejection rates. First, 
in model A we observe that for p = 0.01 the system self- 



The myopic adaptive setting has improved remarkably the 
resilience to congestion without the need of tuning any ex- 
ternal parameters. However, the existence of an abrupt phase 
transition, again as found in I39ll42[l43[l48ll . demands for fur- 
ther improvements. The main goal in order to soften such 
abrupt transition is to avoid that all the nodes close their doors 
due to its own congestion by incorporating an empathetic be- 
havior based on the local knowledge about the dynamical state 
of their neighbors. This empathetic behavior should motivate 
congested nodes to open their doors when detecting an hyper- 
congested state in its surroundings. To this aim we take model 
B 115 511 and reformulate its equations as follows: 

^ +1 = ^ + /3[(l-a)Al- a (A*.) ri ] • (13) 

In the above equations we introduce a new term accounting 
for the average level of congestion in the neighborhood, Tu of 
a node i, 

N . 

The relative importance that nodes assign to the local level of 
congestion in their neighborhoods with respect to their own 
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FIG. 5: (Color online). (Top) Phase diagram p(p) of the empathetic 
routing model for several values of the empathy parameter a. The 
phase diagram of the minimal routing model (circles) is also plot- 
ted for the sake of comparison. (Bottom) For the same empathy pa- 
rameters we show the average rejection rate (77) as a function of p. 
The function (eta) (p) obtained analytically from global minimiza- 
tion and computed from equation l|9} is also shown. 



state is controlled by the parameter a. In particular, when a = 
we recover the myopic setting whereas for a = 1 routers 
behave "altruistically" and their decisions are based solely on 
their neighbor's state of congestion. Thus, the parameter a 
measures the degree of empathy between connected routers. 

In the top panel of Fig. [5] we plot the phase diagrams 
for several values of a together with that of the minimal 
non-adaptive routing model. We observe that for a < 0.5 
the phase-transition is similar to that of the myopic adaptive 
model (a = 0), i.e. showing a critical load of p c ~ 0.1 fol- 
lowed by a first-order transition to full congestion. However, 
from the figure we observe that when a > 0.5 the transition 
to congestion occurs smoothly, thus recovering the behavior 
of the minimal model. On the other hand, the value of p c also 
decreases with a (thus anticipating the onset of congestion) 
although it remains close to the original value p ~ 0.1 until 
a ~ 0.63. Moreover, forp > p c , the curves corresponding to 
a = 0.63 and a = 0.75 reach levels of congestion similar to 
those observed in the minimal model. 



In order to gain more insight about the strategy adopted in 
the empathetic setting we have computed the average level of 
rejection rate as a function of p for the relevant values of a. In 
the bottom panel of Fig. [5]we observe that those curves corre- 
sponding to a > 0.5 are quite different from those obtained in 
Fig. [3] for the myopic adaptive setting. In particular, when 
p « p c the empathetic adaptability shows a large amount 
of rejection. However, as p increases the average rejection 
rate decreases monotonously. This high-rejecting behavior for 
p < Pc, was not observed in the myopic scheme. Quite on 
the contrary, it was shown that nearly all the doors were open 
in the sub-critical regime. However, the high rate of rejec- 
tion observed in Fig. |5]is due to the large degree of empathy 
(a > 0.5) and the existence of a number of under-congested 
nodes, A; < 0, in the sub-critical regime. Under these low 
traffic conditions, most nodes will close partially their doors 
when detecting under-congested neighborhoods, (A*-)r 4 < 0, 
in order to benefit from the availability of neighbors to handle 
their packets. This situation is highly dynamical and most of 
the nodes experiment large fluctuations in their rejection rates 
until the system equilibrium is reached. This microscopic sce- 
nario, although clearly different from that of the myopic set- 
ting, enables to delay the onset of congestion in an efficient 
way. On the other hand, as p approaches p c and for p > p c 
we observe that (for a > 0.5) the value of (77) decreases to 
as p increases. This is due to both the large number of over- 
congested neighborhoods, (A*-)r 4 > 0, surrounding routers 
in the super-critical regime and their large degree of empathy. 
As expected, empathy prevents from the sudden door closing 
when p > p c , thus favoring a smooth phase transition display- 
ing congestion levels simlar to those observed in the minimal 
routing model in the super-critical regime. 

Interestingly, the monotonous decrease of (r?)(p) from 
(77) = 1 at p = shown in Fig. [5J points out a similar be- 
havior to that obtained by means of global minimization of 
congestion. As shown in the bottom panel of Fig. |5]the theo- 
retical estimation of (77) (p) (circles) follows the same trend as 
the self-adopted patterns for a > 0.5. To analyze in detail the 
similarity between the empathetic setting and the microscopic 
patterns predicted by global minimization of congestion we 
plot in Fig. |6]the average value of the rejection rate as a func- 
tion of the degree k of the nodes, (rj)(k), for several values 
of p and a. The panels correspond to (a) p = 0.02 (free- 
flow regime), (b) p = 0.1 (critical point), (c) p = 0.3 and (d) 
p = 0.6 (congested state). The shape of each curve (rj)(k) be- 
haves similarly to the theoretical one as predicted from equa- 
tion (O. More importantly, for each value of p there is one 
value of a, a° pt , for which the curve (77) (fc) fits perfectly the 
prediction made by global minimization of congestion. The 
precise value of a opt depends on p. In particular, forp = 0.02 
we find a opt ~ 0.63, forp = 0.1 we obtain a opt ~ 0.55, for 
p = 0.3 we have a opt ~ 0.68 and, finally, for p = 0.6 the 
value found is a opt ~ 0.75. Moreover, from the top panel of 
Fig. [5] we observe that the values found for a° pt are those 
for which congestion, p(p), is minimum. This result points 
out that empathetic adaptability is able to avoid congestion by 
means of only local information as much as global minimiza- 
tion does. 
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FIG. 6: (Color online). Distribution of the mean rejection rates (77) across degree-classes of the empathetic adaptive model for several values 
of the empathy parameter a compared with the global minimization prediction (dashed line). Different traffic values p are presented: (a) 
p = 0.02 (free-flow regime), (b) p = 0.1 (critical point), (c) p = 0.3 and (d) p = 0.6 (congested state) 



VI. CONCLUSIONS 

We have studied a novel mechanism that allows routers to 
adapt their individual strategies based on their local knowl- 
edge about congestion. Although in our approach nodes can 
only decide either to refuse or to accept incoming packets 
from their first neighbors, we obtain a variety of dynamical 
behaviors. First, we have analyzed the situation when no in- 
dividual adaptability is allowed. This allows us to show that 
whenever a small level of rejection is applied indistinctly to 
all the nodes, one obtains a worse overall behavior than when 
all incoming flows are accepted by the routers. Then, we 
have considered that routers can have different rejection rates 
and derived analytically their patterns to minimize congestion, 
considering global knowledge of the network topology. With 
these globally optimized patterns the resilience to congestion 
of the system can be enhanced significantly. Besides, these 
patterns reveal a dependence of the rejection rate and the de- 
gree of the router while its mean value decays with the incom- 
ing load of packets. 

After deriving global minimization of congestion we have 
studied the situation in which nodes self-adjust their own re- 
jection rates dynamically depending on their instant level of 
congestion (myopic setting). In this case we have shown that 
the critical load of the network is shifted to a value similar 
to that found analytically by means of global minimization of 
congestion. This improvement is again achieved by a proper 
distribution of the rejection rates according to the degrees 
of the routers. However, in the adaptive case, such degree- 
correlated configuration is self-tuned by the system and differs 
from that obtained analytically. As usual in congestion-aware 



routing schemes, such delay in the congestion onset comes 
together with an abrupt transition from the free-flow phase to 
the congested one that prevents from having any warnings of 
the approach to the onset of congestion. For this reason, we 
have finally explored the situation in which routers also con- 
sider the congestion state of their first neighbors to adapt their 
rejection rates. We have shown that when nodes empathize 
with the congestion state of their neighbors, thus not rejecting 
packets from them when they detect an over-congested neigh- 
borhood, the shift in the critical load (obtained through global 
minimization and the myopic adaptability) is preserved and 
followed by a smooth congestion transition. Moreover, the 
analysis of the microscopic patterns of rejection rates when 
empathy is the mechanism at work points out a similar orga- 
nization to that obtained from global minimization. In partic- 
ular, it is possible to find the degree of empathy that perfectly 
agrees with the analytical estimation of the rejection pattern 
that minimize congestion for a given load of information. 

In summary, we have shown that allowing routers to adapt 
their own strategies together with a certain degree of local em- 
pathy is strongly beneficial to the behavior of complex com- 
munication systems. Moreover, the improvement shown when 
local empathy is at work is similar to that obtained by min- 
imizing congestion by means of a global knowledge of the 
network topology. Thus, the empathetic setting represents, a 
remarkable example of how local rules can achieve levels of 
functioning as optimal as those obtained with global knowl- 
edge of the system. Besides, our results open the relevant 
question about how local empathy can be naturally tuned as a 
function of the external inputs. 
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